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ABSTRACT 

Context. Red giant stars are important tracers of stellar populations in the Galaxy and beyond, thus accurate modeling of their structure 
and related observable properties is of great importance. Three-dimensional (3D) hydrodynamical stellar atmosphere models offer a 
new level of realism in the modeling of red giant atmospheres but still need to be established as standard tools. 
Aims. We investigate the character and role of convection in the atmosphere of a prototypical red giant located close to the red giant 
branch (RGB) tip with atmospheric parameters, T eS = 3660 K, logg = 1.0, [M/H] = 0.0. 

Methods. Differential analysis of the atmospheric structures is performed using the 3D hydrodynamical and ID classical atmosphere 
models calculated with the CO 5 BOLD and LHD codes, respectively. All models share identical atmospheric parameters, elemental 
composition, opacities and equation-of-state. 

Results. We find that the atmosphere of this particular red giant consists of two rather distinct regions: the lower atmosphere dominated 
by convective motions and the upper atmosphere dominated by wave activity. Convective motions form a prominent granulation 
pattern with an intensity contrast (~ 18%) which is larger than in the solar models (~ 15%). The upper atmosphere is frequently 
traversed by fast shock waves, with vertical and horizontal velocities of up to Mach ~ 2.5 and ~ 6.0, respectively. The typical diameter 
of the granules amounts to ~ 5 Gm which translates into ~ 400 granules covering the whole stellar surface. The turbulent pressure in 
the giant model contributes up to ~ 35% to the total (i.e., gas plus turbulent) pressure which shows that it cannot be neglected in stellar 
atmosphere and evolutionary modeling. However, there exists no combination of the mixing-length parameter, a M LT, and turbulent 
pressure, P tU rb, that would allow to satisfactorily reproduce the 3D temperature-pressure profile with ID atmosphere models based on 
a standard formulation of mixing-length theory. 
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1. Introduction 

Convection plays an important role in governing the interior 
structure and evolution of red giants (i.e., stars on the red and 
asymptotic giant branches, RGB/AGB). Besides of aiding the 
energy transport from the stellar interior to the outer layers, 
convection mixes heavy elements from the nuclear burning lay- 
ers up into the stellar envelope and atmosphere. Since convec- 
tive mixing changes the local chemical composition it alters the 
stellar structure because of changes in the opacities, thermody- 
namic properties, and nuclear reaction rates of the stellar plasma. 
This affects the observable properties of a star, and a proper 
understanding of convection is thus of fundamental importance 
for building realistic models for stellar structure and evolution, 
which, in turn, are fundamental building blocks of our under- 
standing of individual stars and stellar populations. 

Convection in current one-dimensional (ID) hydrostatic stel- 
lar atmosphere models is treated in a simplified way, typically, 
using the classical mixing-length theory (MLT, Bohm-Vitense 
1958 ) or one of its more advanced variants (Canuto & Mazzitelli 



1991 



Canuto et al.|19 96 ). This approach has a number of draw- 



backs. For instance, the efficiency of convective transport in the 



framework of MLT is scaled by the apriori unknown mixing- 
length parameter, #mlt- It * s commonly taken as a fixed ratio 
of the mixing-length to the local pressure scale height, and is 
usually calibrated using theoretical models of the Sun. Since 
the MLT is a rather simplistic approach, #mlt needs not to be 
the same in main-sequence stars, subgiants, giants and super- 
giants, as is norm ally assumed in the calculation of stellar at- 
mosphere models ( Castelli & Kurucz 2003 ; Brott & Hauschildt 
2005 ; Gustaf sson et al.|2008|) or stellar evolutionary tracks and 
isochrones jD emarqu e et aL] |2004[ jVa ndenber g et al.| [2006 ; 
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Potter et al.|2008HBertelli et al.|2008||2009F 

Despite significant efforts made during the last few decades 
to improve the treatment of convection in stellar structure models 
a fundamental breakthrough is still missing. Seeking numerical 
solutions of the underlying radiation-hydrodynamical equations 
is a promising way to make progress since such models lay-out 
a clear path towards increasing realism in the description of con- 
vection. In this class of models three-dimensional (3D) hydro- 
dynamical model atmospheres have already demonstrated their 
excellent capabilities of reproducing the observed properties o f 
surface convection in the Sun (e.g., Stein & Nordlund 1998), 
and associated spectral d iagnostics (e.g., |Asplund et al.||2000| 
Caffau et al.||2008| |2010|), in many cases outperforming clas- 
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sical ID models despite lacking tunable parameters as in ID 
Similarly successfully, such 3D models wer e applied to other 
types of stars, such as late-type dwarfs ([Gonzalez Hernandez et 



al.|2009[|Ramirez et al.|200 9; Behar a et al.|2010fr and subgiants 
( [Collet et al.|2009| ). 

The progress in the 3D modeling of red giant atmospheres 
has been considerably slower. To large extent this is related to 
the fact that giant models are computationally more demanding, 
and the calculation of a grid of giant atmosphere models is a 
sizable task. Among the early eff orts, a 3D model atmo sphere 
of a red giant was discussed by Kucinska s et al. ( |2005| ). The 
authors have compared broad band photometric colors as pre- 
dicted by the classical ID and 3D hydrodynamical model atmo- 
spheres and have shown that for certain color indices the 3D- 
1D differences may reach -0.25 mag. A detailed analysis of the 
spectral line formation in the atmospheres of somewhat warmer 
red giants (r eff * 4700 - 5100K, logg = 2.2 , [M/H] = 0.0 
to -3.0) was carried out by Collet et al. ([2007). An interesting 
finding of this work is that convective motions may produce sig- 
nificantly cooler average temperatures in the outer atmospheric 
layers, an effect which is increasingly pronounced at low metal- 
licities ([Fe/H] < -2.0). Similar effects have been seen in the 3D 
atmosphere models of late- t ype dwarfs too, see, e.g., Asplund 
|et al,| fll999T ), |Behara et al] ( |2010| ), |Gonzalez Hernandez efaT 
( 2010| ). Spectral lines of various atoms and molecules appear 



typically stronger in a 3D than in a ID model which leads to 
a disagreement between the chemical abundances, e(X), reach- 
ing A log e(X) = -0.5 ... - 1.0 dex at the lowest metallicities. 
Despite these successful attempts, 3D hydrodynamical model at- 
mospheres of red giants still need to be consolidated including 
comp arisons to observations as exemplified by Rami rez et al. 
( |2010| ) who performed a study of a metal-poor giant with em- 
phasis on hydrodynamic properties. 

To further broaden the 3D model basis of giants, we under- 
took a study of the influence of convection on the structure and 
observable properties of red giants. The models were calculated 
with the 3D radiation-hydrodynamics code CO 5 BOLD and will 
eventually cover the entire range of stellar parameters typical for 
stars on the red and asymptotic giant branches (RGB and AGB, 
respectively). Some of these models are already availabl e as part 
of the CIFIST grid of C0 5 B0LD 3D model atmospheres ( |Ludwig| 



et al. 2009). This homogeneous set of 3D model atmospheres 



is well suited to investigate the role of convection in the atmo- 
spheres of red giants of different effective temperatures, gravities 
and metallicities. 

This paper summarizes the first results of the project, fo- 
cusing on the role of convection shaping the atmospheric struc- 
ture of a solar-metallicity red giant located close to the RGB 
tip. The atmospheric parameters of the model are T e s =3660 K, 
logg=1.0, [M/H]=0.0. The analysis is done differentially by 
comparing 3D hydrodynamical and ID static model atmo- 
spheres calculated for the same set of atmospheric parameters 
and with identical opacities and equation-of-state. In two com- 
panion papers we will further discuss the effects of convection 
on the observable properties of this particular red giant, by tak- 
ing a closer look at the formation of individual spectral lines and 
global properties of the spectral energy distribution (Kucinskas 
et al. 2012a,b , in prep.) The present paper mostly describes mor- 
phological properties of the particular model which we, however, 
generalize in places to make statements about convection in red 
giants in general. 



2. The model setup 

Two stellar atmosphere models were used in this study to as- 
sess the influence of convection on the atmospheric structures 
of a red giant: 3D hydrodynamical and classical ID models, 
calculated with the codes CO 5 BOLD and LHD, respectively. Both 
models were computed using identical atmospheric parameters, 
T eff * 3660 K, logg = 1.0, [M/H] = 0.0. According to theo- 



retical evolutionary tracks (e.g., |Cordier et al.||2007| ) this set of 
atmospheric parameters characterizes a late-type giant located 
close to the RGB tip, with a mass of ~ 2 M and age of ~ 1 Gyr. 
In turn, the mass and gravity give a radius of ~ 75 R©. 



2.1. The 3D CO 5 BOLD model 

The 3D atmosphere model was calculated using the CO 5 BOLD 
stellar atmosphere code. The CO 5 BOLD code employs a Riemann 
solver of Roe type to integrate the equations of hydrodynam- 
ics and simultaneously solve for the frequency-dependent radi- 
ation field on a 3D Cartesian grid ([Freytag et al.||2002| |2003 



Wedemeyer et al. 2004 Freytag et al. 2010). For a detailec 



description of the code and its applications see Freytag et al. 
( |20T2l ). 



The code was employed in the 'box-in-a-star' set-up using a 
grid of 150 x 150 x 151 mesh points (x x y x z, where z is verti- 
cal dimension), with a corresponding size of the computational 
domain ofl5.6xl5.6x 8.6 Gm 3 . Effects related to the stellar 
sphericity are neglected. We applied open boundary conditions 
in the vertical direction and periodic boundaries in the horizon- 
tal direction. The convective flux at the lower boundary was con- 
trolled by specifying the entropy of the inflowing gas. The model 
spans a range in optical depth of -6.5 < log <tr 0S s) ^ 5 and sam- 
ples the atmospheric layers over ~11 pressure scale heights. 

In the model calculations we used monochromatic opacities 
from the MARCS stellar atmosphere package ( Gustaf ssonetaL 
2008), grouped in to 5 opacity bins (for more on opacity group- 
ing scheme s ee [N ordlund 1982; Ludwig 19 92} [Ludwig et aL 
|1994[ [Vpgler 2004]). Opacities wer e calculated assuming sola r 
elemental abundances according to Grevesse & Sauval (1998 ), 
with the exception of carbon, nitrogen and oxygen, for which the 
following values were used: A(C)=8.41, A(N)=7. 8, A(Q)=8.67 
simil ar to the CNO abundances recommended by Asplund et al. 
(2005). Local thermodynamical equilibrium (LTE) was assumed 
throughout the entire atmosphere and scattering was treated as 
true absorption. 

The equation-of-state (EOS) used in the model simulations 
takes into account the ionization of hydrogen and helium, as 
well as formation of H2 molecules according to Saha-Boltzmann 
statistics. The ionization of metals is ignored in the current ver- 
sion of EOS since its importance for the gross thermodynamical 
properties is minor. 

After initial relaxation to a quasi- stationary state, the model 
simulations were run to cover a span of ~ 6 x 10 6 sec (~ 70 days) 
in stellar time. This corresponds to ~ 7 convective turnover times 
as measured by the Brunt- Vaisala* and/or advection timescales, 
the latter equal to the time needed by the convective material to 
cross one mixing-length, A = #mlt#p, where H v is the pres- 
sure scale height (both timescales were estimated at tr oss — 1.0, 
see Sect. |3.2.3| ). A set of 70 relaxed 3D model snapshots was 
used in our analysis of the atmospheric structure. This is suffi- 
cient for the aspects addressed in this work. However, one should 
keep in mind that the overall statistics gathered is limited (in part 
dictated by computational cost), and the precision of convection- 
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Fig. 1. Snapshots of the emergent white light intensity during the temporal evolution (ordered upper-left to lower-right) of the 
hydrodynamical red giant model atmosphere. Bright granules where the matter is rising are surrounded by darker and significantly 
narrower intergranular lanes associated with down-flows. The spatial size of each frame is 15.6x15.6 Gm 2 , and the time interval 
between frames 300 ks. The relative RMS intensity contrast is given in the upper left corner of each frame. 



related properties (e.g., the turbulent pressure) is likely not better 
than -5%. 



2.2. The average (3D) C0 5 B0LD model 

Since the 3D hydrodynamical and ID classical model atmo- 
spheres are based on different physical assumptions, one may 
naturally expect differences in their resulting properties. Two as- 
pects are worth mentioning here. First, the 3D hydrodynamical 
model atmospheres predict significant horizontal fluctuations in 
temperature, density, pressure, velocity and other kinematic and 
thermodynamic properties. This may produce horizontal varia- 
tions in, e.g., shapes and strengths of spectral lines formed in 
different parts of the model photosphere. Second, convection is 
a natural process arising in hydrodynamical model atmospheres 
inherent to the equations of hydrodynamics and radiative trans- 
fer. Besides convection as such, hydrodynamical models ex- 
hibit significant overshoot of material into the upper layers of 
the atmosphere which should be convectively stable under the 
Schwarzschild criterium. None of these effects is properly ac- 
counted for in the ID classical models built on the prescription 
given by MLT. The consequence is that the structure of 3D hy- 
drodynamical models, even if horizontally averaged, is different 
from ID models calculated with identical atmospheric parame- 
ters. 



The relative importance of these effects can be assessed 
using the average 3D model. Such a model can be calcu- 
lated by horizontally averaging the thermal structure over sev- 
eral instances in time (snapshots). By definition, the new one- 
dimensional structure, further referred to as (3D) model, retains 
no information about the horizontal inhomogeneities but keeps 
the imprints from the different (i.e., hydrodynamical) treatment 
of convection in the original 3D model. The 3D-(3D) and 
(3D)- ID differences may therefore provide valuable informa- 
tion about the relative importance of horizontal fluctuations and 
the time-dependent effects related with the treatment of convec- 
tion, respectively. 

The (3D) model used in this work was obtained by averag- 
ing 70 3D snapshots over the surfaces of equal Rosseland op- 
tical depth. To preserve the radiative properties of the original 
3D model we averaged the fourth moment of temperature and 
first moment of gas pressure, following the prescription given in 
ISteffen et al.| ( [T995] ). 

2.3. 1D LHD model 

Comparison ID hydrostatic models were calculated with the 
LHD code, using the same atmospheric parameters, elemental 
abundances, opacities and EOS as in the 3D model calcula- 
tions described above. Convection in the LHD models is de- 
scribed using the mixing-length theory, adopting the formulation 
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Fig. 2. Non-normalized probability density of the emergent 
white light intensity extracted from 14 snapshots of the 3D 
model, "all" labels the overall distribution, "up" and "down" are 
the distributions restricted to up- or down-flowing material only. 
Intensities are given as relative deviation from the (temporal and 
spatial) mean intensity. As velocity criterion the sign of the ver- 
tical velocity component at log T ross = was taken. 



of Mihalas (1978]). The LHD models characterized by several dif- 
ferent mixing-length parameters were used in this work, with 
(#mlt = 1.0, 1.5 and 2.0. It should be stressed though that the 
choice of the mixing-length parameter affects mostly the deeper 
atmospheric layers and has limited influence on the formation of 
the em itted spectrum in this particular red giant ( Kucinskas et al.| 
|2005] Kucinskas et al. 2012, submitted to A&A). 

In a few cases we added an ad-hoc overshooting in the LHD 
models. This was done by preventing the convective velocity go- 
ing to zero in the formally stable regions but forcing it to a pre- 
scribed, fixed value. The convective flux was calculated from 
the standard MLT formulae. In the overshoot region this results 
in a downward directed convective flux, and usually a flux di- 
vergence which leads to additional cooling. This is intended to 
mimic the cooling effect by overshooting often observed in 3D 
models in the upper photosphere. 

3. Results and discussion 



3.1. Properties of surface granulation 

The 3D model of red giant predicts the existence of surface gran- 
ulation which is clearly visible in the time series of emergent 
white light intensities shown in Fig. [T] Although similar gran- 
ulation patterns a re routinely seen in the 3D mod els of the Sun 
dNordlund|1982|), M- dwarfs ( [Ludwig et al.|2002|), white dwarfs 
( Freytag et al.||1996| ) , brown dwarfs ([Freytag et al.||2010|), pre- 
mai n seq uence stars ( [Ludwig et al.|2006|), sub giants ; ( [Ludwig et| 
|al.|20 06) and warm giants ( Collet et al.|2007| ), the existence of 
granulation in cool red giants is not exactly self-evident: MLT 
predicts that the surface convective zone is confined to opti- 
cally thick layers in the ID models of giants, with an only thin, 
marginally unstable region extending into the upper atmosphere 
(see the ID stratification of the convective velocity in Fig.fTO]). 




Fig. 3. Emerging white light intensity (top panels), vertical (v z ) 
and horizontal (vh = ^vl + Vy) velocities (bottom panels) in a 

typical 3D snapshot of the red giant model, all at log tr oss = 0. 
Yellow (bright) and brown (dark) pixels in the top-right and bot- 
tom panels highlight the areas where the relative intensity devia- 
tion is larger than AI/(I) = 0.3 in the up-flows and down-flows, 
respectively. Locations where the vertical velocity v z < -1 Mach 
(supersonic down-flows) are marked with white contours, those 
where the vertical and horizontal velocities exceed 1 Mach are 
contoured in black (contours start at Mach = +1.0 and continue 
at +0.5 Mach thereof). The relative white light intensity contrast 
is 18.2 % for this particular snapshot. 



boundary of the convective region and overshoots into the upper 
atmospheric layers. Similarly to dwarfs, the efficiency of convec- 
tive overshoot exhibits an exponential decline with height, e.g., 
as hinted at by the shape of the vertical velocity profile in the 
op tical depth range log t r oss ~ 1.0 ... - 2.0 (Fig. 10). According 
Freytag et al.| ( |1996| ), such exponential decline is caused by 



Consequently, the outer optically thin regions are essentially un- 
affected by convection in ID models. 

In the 3D model, the geometric distance between the upper 
boundary of the convective region and the optically thin region 
is, however, not large. The material, therefore, crosses the formal 



to_ 

the penetration of convective modes with long horizontal wave- 
lengths into the formally convectively stable layers. 

It may perhaps be interesting to note that intergranular lanes 
are narrower in our giant model than those seen in the 3D models 
of the Sun but wider than in the models of late-type dwarfs (cf. 
|Ludwig et al.||2006] ). Why this is so is not immediately clear, 
although the increase of the relative width of down-drafts with 
T e ff suggests that it may be related to a stronger smoothing of 
thermal inhomogeneities caused by the more intense radiative 
energy exchange at higher temperatures. 

Despite the fact that the granulation pattern shown in Fig. [T] 
exhibits the typical appearance known from the 3D model atmo- 
spheres of other types of stars, Fig. [2] illustrates that its intensity 
distribution does not show the familiar bimodal shape related to 
the dark and bright areas of the granulation pattern but a sin- 
gle maximum. The intensity distributions restricted to up- and 
down-flowing regions makes it clear that the correlation between 
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hot up-flowing and cool down-flowing gas is still present, a clear 
indication that convection as such takes place. 

However, the intensity distributions in Fig. [2] illustrate that 
this correlation is far from perfect. For instance, high surface 
brightness does not necessarily imply that matter is in the up- 
flow since some of the down-flows may also be very bright. As 
seen in Fig. [3] such bright down-flows are typically seen on the 
edges of convective cells and in many cases are located imme- 
diately next to the uprising material that is also characterized by 
high white light intensities. The high intensity in these regions 
is caused by their significantly higher temperatures (see Fig. [5}. 
These are, in turn, caused by the dissipation of fast horizontal 
flows (mostly weak shock waves) when they collide with the 
down-flowing material at the granule edges are and are deflected 
downwards. This also explains why in some cases (but not al- 
ways) the brightest down-flows are located immediately next to 
the brightest up-flows. There are also isolated regions of bright- 
est intensity that are located within the granules and which are 
tracing the hottest uprising material. 

Similar edge-brightened granules are seen on the surface of 
the Su n too, both in observ ations (e.g.,|Keller & von der Luehe 
T992] ) and 3D models (e.g., |Stein & Nordlund|1998| ). However^ 
horizontal and vertical velocities are significantly higher in the 
atmosphere of the red giant model. While only mild shocks are 
seen in the models of the Sun, with the maximum Mach num- 
bers of ~ 1.5 and ~ 1.8 in the vertical and horizontal directions, 
respectively, the corresponding numbers in the outer atmosphere 
of red giant may reach to ~ 2.5 and ~ 6.0 (Fig. [4]). It should also 
be mentioned that some of the brightest spots in the up-flows 
sometimes appear not on the edges but closer to granule cen- 
ters, delineating the regions that are splitting into new granules 
(Fig.|3J. 

The average intensity contrast of the granular pattern seen 
in our model red giant is -18.1% (calculated for a sequence of 
70 3D snapshots). This is somewhat higher than the intensity 
contrast in the 3D atmosphere models of the Sun (~ 15%) but 
slightly lower than the corresponding number in a model of a 
subgiant (~ 23%, |Ludwig et al.|2002| ). 

The size of a typical granule in the red giant model is of the 
order of ~ 5 Gm (Fig.[T]), somewhat larger than 10 times the pres- 
sure scale height at the surface - depending on the exact location 
and whether turbulent pressure is considered in its definition or 
not. The relative size is on the high side in comparison to mod- 
els of higher gravity. The high horizontal velocities found in the 
model are in part a consequence of this large granular size. In 
absolute terms, the size is by at least three orders of magnitude 
larger than the typical size of sol ar granules, roughly equal to 
~ 1 Mm (e.g., Nordlund et al. 2009| ). The latter number translates 
into ~ 2 x 10° granules on the Sun, in contrast to only -400 
granules on the surface of the red giant studied here. This, to- 
gether with the high intensity contrast, indicates that granulation 
causes larger temporal fluctuations of the observable properties 
of giants than of dwarfs, in the simplest case of their brightness. 



3.2. Thermal structure of the atmosphere 
3.2.1 . Properties of the full 3D structures 



As it was already discussed in Sect. 3.1 the 3D model shows a 
prominent granulation pattern, a direct consequence of convec- 
tive motions. Convective cells are discernible in the deeper (and 
hotter) atmospheric layers and are pronounced between the op- 
tical surface (at z = cm or <tr oss ) = 1) and the lower boundary 
of the model at z = 7 x 10 11 cm (Figs. and |5}. In this depth 



range convection manifests itself in the form of wide up-flows 
and narrower and cooler down-flows. 

Velocity amplitudes increase in the higher atmospheric lay- 
ers and their vertical velocity, v z , sometimes becomes supersonic 
(upper panels in the left column of Fig. [4]). The granulation pat- 
tern looses coherence towards the upper atmosphere, and the 
velocities become gradually dominated by motions related to 
acoustic waves. 

The horizontal flow speed, Vh, is low in the subphotospheric 
layers due to the predominantly vertical motions here. The situ- 
ations starts to change when matter approaches the optical sur- 
face. The decreasing opacity leads to enhanced photon losses, 
the matter rapidly cools, becomes denser, and its vertical ve- 
locity decreases. A density inversion is formed just below the 
optical surface, at log tr oss « 0.5 (see Fig. [9]). Simultaneously, 
the up-flowing material is deflected sideways until it reaches the 
granule edges finally merging into the down-flows. The result- 
ing pattern of horizontal velocities is tracing the granular shapes 
and is visible in the atmospheric layers above the optical surface 
(Fig.gJ). 

Figure [6]illustrates that the amplitude of the velocity fluctua- 
tions on horizontal planes is highest in the outermost layers, and 
largely shaped by shock waves. This holds for both the vertical 
and horizontal velocities. Towards deeper layers, the fluctuations 
in the vertical velocity become noticeably smaller, fluctuations 
of the horizontal velocity also decrease but to a lesser degree. 
As alluded to before, this is a consequence of the take-over of 
the more regular convective motions over wave motions. After 
passing a minimum around the optical surface the fluctuations 
of the vertical velocity increase again, a signature of the roughly 
columnar flow pattern in the the subphotospheric layers. 

As shown in Fig. [7] the columnar pattern in the subphoto- 
sphere is also imprinted in the temperature field which is rather 
homogeneous in the up-flows, and drops in the temperature 
marking down-flows. The temperature fluctuations reach their 
maximum in the surface layers, and are quickly reduced in am- 
plitude in the optical thin layers. This is surprising in view of the 
substantial density fluctuations (see Fig. |7J lower panel) present 
in the photospheric layers. This indicates an efficient smooth- 
ing of the temperature field by radiative energy exchange coun- 
teracting temperature changes by adiabatic compression or ex- 
pansion. The vertical cuts in Fig. [8] and [9] show that the shock 
waves form a rather irregular pattern in the upper atmosphere. 
The shock fronts are often horizontal or arc-like, similar to 
those seen in the 3D hydrodynamical models of the Sun (e.g., 
|Wedemeyer et al.|2004b . Most of the shock activity takes place 
above z = -0.5 x 10 a1 cm (logTR OSS < -1.0) where the flow 
density is low enough for the flow to accelerate to supersonic 
speeds. 

Finally, we would like to point out that our model is rather 
shallow in the sense that its extent below the optical sur- 
face supercedes the horizontal size of granular cells only little. 
Consequently, the merging and narrowing of downdrafts with 
depth as discussed by Stein & Nordlund (1998) is not clearly dis- 
cernable. However, qualitative similarity of the convective mor- 
phology makes us belief that models with larger extent in depth 
will exhibit the same feature. 



3.2.2. Properties of the (3D) stratification 



As we have already seen in Sect. |3.2.1| two regions can be 
distinguished in the photosphere of the giant model: the lower 
(and hotter) part dominated by convective motions per se and 
upper layers where the flow is driven by acoustic waves. This 
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Fig. 4. Vertical (left) and horizontal (right) velocity maps in a Fig. 5. Temperature (left) and density (right) maps in the typi- 

typical 3D snapshot of the red giant model constructed at dif- cal 3D snapshot of the red giant model constructed at different 

ferent geometrical depth, z, with z - set at log tr oss = and geometrical depth, z, with z = set at log tr oss = and increas- 

increasing towards the stellar center. Locations where the verti- ing towards the stellar center. Contours are drawn at every 200 K 

cal velocity v z < -1 Mach (supersonic down-flows) are marked (panels 1-3, top-down) and 500 K (panels 4-6) for temperature, 

with white contours, those where v z and Vh exceed 1 Mach are and 0.5 dex (panels 1-4, top-down) and 0.1 dex (panels 5-6) for 

lontoured in black (contours start at Mach = +1.0 and are drawn density. Vertical lines mark the x position (x = 7.77 x 10 11 cm) at 

at every ±0.5 Mach thereof). Vertical lines mark the x position which the 1-dimensional velocity profiles shown in Fig. [9] were 

(x = 7.77 x 10 11 cm) at which the 1-dimensional velocity pro- constructed, 
files shown in Fig. [8] were taken. 



H.-G. Ludwig and A. Kucinskas: 3D hydrodynamical CO 5 BOLD model atmospheres of red giant stars. I. 





5 









-5 




5 









-5 




5 







1 


-5 




5 









-5 




5 









.5 




5 









.5 




10 




5 









10 




5 









10 


ST 


5 


E 







10 




5 









10 




5 









10 




5 










0.4 0.8 

x [x 10 12 cm] 



1.6 



Fig. 6. 1 -dimensional vertical (top) and horizontal (bottom) ve- 
locity profiles in the red giant model at x = 7.77 x 10 11 cm 
(marked the by white vertical lines in Fig. [4]) and shown at differ- 
ent geometrical depths as indicated on the left-hand side of each 
velocity profile. Grey horizontal lines indicate the level where 
the velocity is equal to zero. 



substructure clearly manifests itself in the (temporally and hori- 
zontally averaged) RMS vertical and horizontal velocity profiles 
(Fig. [TO]). Convective motions gradually cease beyond the formal 
convective boundary at log tr oss -0.5. However, the matter par- 
tially penetrates beyond it in a form of exponentially decreasing 
overshoot. Qualitati vely, this is very similar t o the picture seen 
in late-type dwarfs ( [Ludwig et al.|2002|[2006| ). 

The temperature profile of the (3D) model follows closely 
the ridgeline corresponding to the maxima of the probability 
density distribution of temperature and optical depth in the 3D 
hydrodynamical model (Fig. [TT]). A sharp decrease both in the 
3D and (3D) temperature profiles between the optical depths of 
log tross = and 2.0 is caused by the rapid increase in the radia- 
tive cooling rate close to the optical surface. 

Interestingly, the temperature profile corresponding to the 
(3D) model is markedly different from the ID temperature pro- 
files in the deeper atmosphere, at logTR OSS > 0.5 (Fig. [T2| ). 
Moreover, none of the ID LHD models with different mixing- 
length parameters, o^mlt, is able to satisfactorily reproduce the 
stratification of the (3D) model. The resulting different pressure- 
temperature relation have a direct influence on the radius of the 
star. Note, that turbulent pressure was neglected in the ID LHD 
models at this stage. 



Figure 13 shows that in the higher atmosphere the 3D model 
exhibits a temperature increase of about 20 K relative to pure 
radiative equilibrium conditions. Energy fluxes show that this 
is not the consequence of mechanical heating (e.g. by shocks). 
Only the radiative flux provides significant energy exchange in 
these layers. LHD models with an ad-hoc overshooting velocity 
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Fig. 7. 1 -dimensional profiles of temperature (top) and density 
(bottom) at x - 1.11 x 10 11 cm (marked by the white vertical 
lines in Fig. [5} and shown at different geometrical depths. 



show a certain degree of cooling. The cooling is not great de- 
spite the substantial overshooting velocities put into the ID mod- 
els. This again illustrates the rather tight coupling of the tem- 
perature to radiative equilibrium conditions as already seen by 
the rather small horizontal temperature fluctuations in the upper 
photosphere. We argue that the heating in the 3D case is due 
to an altered radiative equilibrium temperature in the presence 
of horizontal T-inhomogeneities and the specific wavelength- 
dependence of the opacity (see Appendix [Bj. Interestingly, this 
is opposite to what is seen in the outer atmosphere of red gi- 
ants at lower metallicities, where the average temperature of the 
hydrodynamical model is significantly lower than that of the cor- 
responding ID model (e.g., [Collet et al.|2007 Dobrovolskas et 



|al.|2010t[Fanauskas et al.|2010| ). 



3.2.3. Radiative, hydrodynamical and rotational time scales 

In order to asses the possible interplay between various physical 
phenomena that take place in the red giant atmospheres we cal- 
culated a number of characteristic radiative and hydrodynamic 
time scales (see Appendix |A| for the definitions). The time scales 
were calculated using the 3D model as background, and are plot- 



ted versus the optical depth in Fig. [14 

Two convection-related time scales, as given by the Brunt- 
Vaisala period, £bv, and the time for crossing one mixing-length, 
^adv, are similar within an order of magnitude and show little vari- 
ation throughout the entire atmosphere. The radiative time scale 
taking into account the wavelength dependence of the opacity, 
t m & is significantly larger than convective times scales in the deep 
atmosphere where the flow is nearly adiabatic. However, £ rad 
rapidly decreases with increasing height and above log tr oss ~ 2 
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Fig. 8. 2D profiles of vertical (top) and horizontal (bottom) ve- 
locity in the red giant model at x = 7.77 x 10 11 cm (marked by 
white vertical lines in Fig. [4]). Locations where the vertical veloc- 
ity v z < -1 Mach (supersonic down-flows, top panel) are marked 
with white contours, those where v z and Vh exceed 1 Mach are 
contoured in black (contours start at Mach = +1.0 and are drawn 
at every ±0.5 Mach thereof). 



becomes about two orders of magnitude shorter than either ^v 
or r a d v , thus causing a rapid evolution of the flow towards radia- 
tive equilibrium. For further comparison we added the radiative 
time scale £ ra d(Ross) based on grey Rosseland opacities. In the 
optically thin region it substantially overestimates the radiative 
relaxation time. 

The time scale of thermal relaxation, or Kelvin-Helmholtz 
time scale, t^n is considerably larger than the convective time 
scale deep in the atmosphere (logTR OSS > 4). This may suggest 
that the flow in the deeper atmosphere will not be able to relax 
to thermal equilibrium during the model simulation time, ~ 6 x 
10 6 sec (see Sect.|2]), which is ~ 1-2 orders of magnitude smaller 
than ^kh- However, such conclusion would be wrong, since at 
these depths thermal relaxation is in fact governed by the mass 
exchange due to convection which takes place on significantly 
shorter time scales. 

The assessment of whether stellar rotation may influence the 
dynamical and radiative processes in the atmospheres of red gi- 
ants is not a trivial task, especially since their rotational periods 
are largely unknown. In order to obtain a qualitative picture one 
may nevertheless compare the relevant time scales. To derive an 
estimate of the rotational time scale we resorted to the available 
measurements of projected rotational velocities, v rot sin i. For red 
giants these are are typically in the ran ge of 2 - 9 km/s (e.g., 
|Carney et alT1|2Q08| |Cortes et al.]|2009l ). I snoring the fact that 
rotational velocity, v rot , may depend on the effective tem perature 
( [Carney et al.|2008| ) and metallicity ( [Cortes et al.|2009| ), and tak- 



0.8 

x [x10 12 cm] 



Fig. 9. 2D profiles of temperature (top) and density (bottom) in 
the red giant model at x = 7.77 x 10 11 cm (marked by white 
vertical lines in Fig. [5]). Contours are drawn at every 250 K (top) 
and 0.2 dex (bottom). 
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Fig. 10. Temporally and horizontally averaged RMS vertical ve- 
locity (thick solid line) and RMS horizontal velocity (thick dot- 
dashed line) of the 3D hydrodynamical model in comparison to 
the ID LHD models (thin lines) as a function of Rosseland optical 
depth. The LHD models were calculated using different mixing- 
length parameter, c*mlt = 1-0, 1.5, 2.0 (top-down). 



ing 3 km/s as a representative rotational velocity, together with 
a radius of ~ 75 R© (see Sect. [2]) one obtains the rotation pe- 
riod of ~ 10 8 secQ This is by at least 1-2 orders of magnitude 
larger than the convective and radiative time scales, except in 



1 This compares well with , e.g., the rotational pe riod of 2 + 0.2 yr 
(~ 6.3 x 10 7 sec) derived by Gray & Brown (2006 ) for Arcturus, the 
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Fig. 11. Non-normalized joint probability density function of Fig. 13. Deviations of temperature profiles from the LHD 
temperature and optical depth (thin contour lines), and average (afMLT=2.0) convective-radiative equilibrium model. The green 
temperature profile (thick solid line) of the 3D model. The statis- line shows the (3D) model, red and black lines depict LHD mod- 
tics was obtained from three snapshots of the flow field. The els with different overshooting velocities, v ver=lkms _1 and 
contour lines are plotted in steps of factors of VlO. v ove r=2 km s" 1 , respectively. 




E 



Fig. 12. Temperature profile of the (3D) model (thick solid line) 
in comparison to ID LHD models (thin solid lines) as a func- 
tion of Rosseland optical depth. The LHD models were calculated 
with different mixing-length parameters, c*mlt = 1.0,1.5,2.0 
(upper— slower). 



the deepest layers where the radiative time scale is longer. While 
we therefore conclude that rotation should have a minor influ- 
ence on governing the atmospheric dynamics of this particular 
red giant, it should be noted that in the deeper atmosphere rota- 
tional and convective time scales may become comparable and 
thus the interaction between convection and rotation interaction 
may become non-negligible. 

3.2.4. The role of turbulent pressure 

One of fundamental differences between the 3D hydrodynamical 
and classical ID atmosphere models of red giants is that hydro- 
dynamical models predict non-zero turbulent pressure, P tm b- The 
contribution of the turbulent pressure to the total pressure (i.e., 
gas and turbulent) may be as large as 15% in the 3D hydrody- 
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Fig. 14. Radiative and hydrodynamical time scales in the 3D at- 
mosphere model of red giant. 



namical models of the Sun ( Houdek 2010 ). In order to assess its 
importance in the red giant atmosphere we calculated the turbu- 
lent pressure as P t urb = (p^), where p and v z are gas density and 
velocity, respectively, and angular brackets denote temporal av- 
eraging as well as averaging on surfaces of equal optical depth. 
The ratio of turbulent to total pressure is plotted versus optical 



depth in Fig. 16 



atmospheric parameters of which are not too different from the red giant 
studied here. 



Although rarely done in the routine calculations of, e.g., 
ID stellar model atmosphere grids, in principle turbulent pres- 
sure can be included in the calculation of classical ID atmo- 
sphere models, in addition to gas and electron pressure. The 
result of such an exercise is shown in Fig. 15 where we plot 
two LHD models calculated with non-zero turbulent pressure 
(^turb - fpv 2 > nere v is convective velocity as given in the frame- 
work of MLT, and / a dimensionless factor, usually f < 1). 
Three additional models shown there were constructed using dif- 
ferent mixing-l ength parameter s, q^mlt = 1.0, 1.5,2.0, utilizing 
a formalism of Mihalas (1978 ) and assuming a vanishing tur- 
bulent pressure, P tur b = 0. As it is evident from Fig. [15] there 



is no reasonable combination of #mlt and / allowing to con- 
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Fig. 15. Temperature stratification in the 3D hydrodynamical 
model of a red giant, as a function of gas pressure, P gas . The 
thick solid line is the profile of the (3D) model, the thin lines 
are ID LHD models calculated using different mixing-length 
parameter ((*mlt = 1.0, 1.5,2.0, from left to right; in all cases 
Ptmb = 0). The two ID models with non-zero turbulent pressure 
are shown by dashed (c*mlt = 2.0 and / = 1.0) and dashed- 
dotted (q^mlt = 2.0 and / = 2.0) lines. 



Fig. 17. Entropy profiles as a function of optical depth. The 
dashed line indicates the value of the asymptotically reached en- 
tropy in the 3D model. Thick solid line is entropy stratification of 
the temporally and spatially averaged (3D) model, thin lines are 
ID LHD models calculated with different mixing-length param- 
eters, #mlt = 1.0, 1.5,2.0 (left-to-right). The two models with 
non-zero turbulent pressure are shown by dashed (c^mlt = 2.0 
and / = 1.0) and dash-dotted (#mlt = 2.0 and / = 2.0) lines. 
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Fig. 16. Turbulent pressure to total pressure versus optical depth 
in the 3D model. 



struct a classical ID model which would reproduce the pressure- 
temperature relation of 3D hydrodynamical model. 

How important then is the contribution of turbulent pres- 
sure in the atmospheres of red giants in general? In the 3D 
hydrodynamical model atmosphere of a red giant studied here 
Aurb provides a non- vanishing contribution to the total pressure 
nearly throughout the entire atmosphere, except for the deepest 
layers where its influence is minor (Fig. 16). The distribution 
of ^turb/CPgas + Aurb) is double-peaked: the first maximum at 
log tross ~ 4 is due to increasing relative importance of connec- 
tive flow, whereas the monotonous rise beyond log tr oss 2 is 

due to wave activity in the upper atmosphere. The peak value s 
are significantly larger than those in the Sun (cf. Houdek|2010 ): 
the contribution of turbulent pressure may amount to ~ 25% and 
-13% due to convective motions close to the optical surface in 
the red giant and the Sun, respectively. 



Obviously, the importance of turbulent pressure can not be 
neglected in realistic modeling of stellar atmospheres and stel- 
lar evolution. One may argue that in the case of the red giant 
studied here the discrepancy between the pressure-temperature 
relations of the 3D hydrodynamical and ID models occurs deep 
in the atmosphere (logr Ross ~ 2) and thus should have a mi- 
nor influence on the formation of emitted spectrum. However, 
this may be of importance for giants at lower metallicities where 
convection can reach into the layers beyond the optical surface. 
Additionally, larger pressure at any given geometrical depth 
would lift the higher-lying stellar layers outwards and may even- 
tually change the stellar radius and luminosity, with direct con- 
sequences on the evolutionary tracks and isochrones. Clearly, is- 
sues related with the proper treatment of turbulent pressure in the 
current stellar evolution and atmosphere modeling of red giants 
warrant further study. 

3.2.5. Effective mixing-length parameter for the 1 D model 
atmospheres 

The potential of the hydrodynamical model atmospheres for the 
cal ibration of the mixing-length parameter, c*mlt, was explored 
by |Ludwig et al. ([1999 ). The authors have shown that the en- 
tropy profiles in the 2D models of solar-type stars posses a 
minimum in the convectively unstable region. While the shape 
and location of this entropy minimum is somewhat different in 
individual spatially resolved profiles, the common property is 
that the entropy gradient becomes very small in the deep at- 
mospheric layers where the adiabatic up-flows dominate. Such 
asymptotic behavior is also characteristic to ID model atmo- 
spheres but the asymptotic entropy values reached deep in the 
atmosphere are different for different mixing-length parameters, 
^mlt- Therefore, the height of the asymptotic entropy plateau in 
the 2D hydrodynamical models can be used to calibrate the #mlt 
for use with ID models ( [Ludwig et al.|1 999 ). Further steps in this 
direction were made by |Ludwig et al.| ( |2002| ) who demonstrated 
that this approach works well also with the 3D hydrodynamical 
model of M-type dwarf. 
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In order to check the feasibility of such approach with the 
red giant model studied here, we utilized the set of LHD mod- 
els calculated in the previous section using different values of 
(Imlt and plotted their entropy profiles versus the optical depth 
in Fig. [T7] Also shown there are the temporally and spatially 
averaged specific entropy profiles of the (3D) model. Clearly, 
the 3D entropy profile shows hints of an asymptotic flattening, 
similarly to what is seen in the ID models. We emphasize that 
horizontally resolved entropy profiles show a plateau at a value 
indicated by s Qm in Fig. 17 Interpolation between the asymptotic 
entropy values found in the ID models with different values of 
the mixing-length parameter and zero turbulent pressure yields 
^mlt ~ 1.8. 

One remark of caution is here again related to the turbulent 
pressure. Evidently, one may get significantly different #mlt val- 
ues for the ID models characterized by different turbulent pres- 
sure, Pturb (Fig- [T7] ). As we have seen in the previous section, 
there is no straightforward procedure to calibrate turbulent pres- 
sure with the aid of (3D) models. Moreover, the net effect of the 
individual spatially resolved entropy profiles will be always dif- 
ferent from that given by the (3D) model. These facts once again 
stress the importance and difficulties of the proper incorporation 
of turbulent pressure into the atmospheric and evolutionary mod- 
els of red giants. 

3.2.6. How realistic is our 3D hydrodynamical model? 

Obviously, the realism of the 3D hydrodynamical model dis- 
cussed in this study is limited by a number of the assumptions 
used in the modeling procedure. Some issues need to be men- 
tioned here. 

Our model of a red giant has limited grid resolution of 
150 x 150 horizontal grid points spread over ~ 1/6 x 1/6 of the 
stellar surface. While such coverage and resolution may be suf- 
ficient to pin down the general properties of atmospheric stratifi- 
cations, substantially better resolution would be needed to inves- 
tigate the properties of atmospheric motions on smaller scales. 
Moreover, the fraction of the surface area covered by the model 
is already large and so the effects of sphericity may become im- 
portant. These issues should be addressed in future work. 

As was described in Sect. 2.1 we used monochromatic opac- 
ities grouped into 5 opacity bins to reduce the computational 
work load. Previous studies and tests that we did in the course of 
the present study have demonstrated that opacity grouping yields 
a surprisingly good agreement with the exact radiative transfer 
calculations obtained using opacity distribution functions, ODFs 
(see, e.g., |Nordlund|1982||Collet et al.|2007] ). One has to be cau- 
tious, however, since the effects related to molecular opacities 
may become increasingly important in the outer atmospheres of 
red giants, and the opacity grouping into small number of bins 
may be insufficient. 

One more issue of concern may be related with the treatment 
of scattering. It has been shown by |Collet et aL] p007| |2008| ) 
that scattering treated as true absorption may produce signifi- 
cantly warmer temperature profiles in the outer atmospheres of 
red giants at low metallicities, compared to those that are ob- 
tained when scattering is treated properly. The treatment of co- 
herent and isotropic scattering in the solution of radiative trans- 
fer equation was recently implemented in the 3D hydrodynam- 
ical code BIFR0ST by |Hayek efaL] ( |20 1 0| ) . While it seems that 
in general the continuum scattering does not have large effect 
on the thermal structure in the atmospheric model of the Sun, 
proper inclusion of scattering in the line blanketing may reduce 
the temperature by ~ 350 K in the optically thin layers below 



logT5ooo ^ _ 4 (Hayek et al. 2010). Obviously, this issue has 
to be addressed properly in the future 3D hydrodynamical mod- 
eling of red giants, however, is mainly an issue for metal-poor 
objects. 

4. Conclusions and outlook 

In this study we have utilized a 3D hydrodynamical model of a 
red giant calculated with the CO 5 BOLD code to investigate the in- 
fluence of convection on the properties of its atmospheric struc- 
tures. The model had solar metallicity and its atmospheric pa- 
rameters, r e ff = 3660 K and log g = 1.0, were typical to those of 
a red giant located close to the RGB tip. We also used a number 
of ID atmosphere models calculated with the LHD code to com- 
pare the predictions of the 3D hydrodynamical model with those 
of classical ID models. The ID LHD model atmospheres shared 
the same opacities and equation-of-state as used in the calcula- 
tions of the 3D hydrodynamical models and covered the same 
range in the optical depth. Thus, the comparison was done in a 
strictly differential way, to ensure that the differences revealed 
would reflect solely the differences related to dimensionality be- 
tween the two types of models. 

Similarly to what is seen in dwarfs, the giant model predicts 
the existence of pronounced surface granulation pattern with a 
white light intensity contrast larger than in solar models (~ 18% 
versus ~ 15%, respectively). The size of the granules relative 
to the available stellar surface is significantly larger in the giant 
than in the Sun, with a typical ratio of ~ 400 and ~ 2 x 10 6 gran- 
ules over surface area, respectively. This, together with the larger 
intensity (and temperature) contrast in the giant model suggests 
that convection may have substantially larger influence on the 
observable properties (such as spectral energy distribution, pho- 
tometric colors) of red giants than in the Sun. Exploratory work 
in this direction seems to confirm this prediction ( Kucinskas et 
|al.|2009l[20T0| ). 

It is important to note that the atmosphere model displays 
a significant activity of shock waves, especially in the outer at- 
mosphere. In comparison with the mild shocks in models of the 
Sun, with vertical and horizontal velocities of up to Mach -1.5 
and -1.8, respectively, the shocks in the atmosphere of the giant 
are considerably stronger, with corresponding Mach numbers of 
-2.5 and ~ 6. We also find that the mean temperature of the hy- 
drodynamical model is slightly higher than that of the ID model 
with identical atmospheric parameters. 

Turbulent pressure is considerably more important in the gi- 
ant model than it is in the Sun and may amount to ~ 35% of 
the total pressure (i.e., gas plus turbulent). This leads to a lifting 
of the outer atmospheric layers to larger radii, besides, it would 
alter the pressure-temperature relation. Interestingly, no combi- 
nation of the mixing-length parameter and turbulent pressure 
in the ID models satisfactorily reproduce the mean pressure- 
temperature relation of the hydrodynamical model. 

We used the 3D hydrodynamical model to obtain the effec- 
tive mixing-length parameter for the classical ID stellar atmo- 
sphere work. Interpolation between the entropy profiles of the 
ID models calculated with different mixing-length parameters 
and zero turbulent pressure yields <*mlt = 1.8. This solution, 
however, is not unique since identical mixing-length parameter 
may be obtained for models characterized by different amount 
of turbulent pressure. The question of how to treat the turbulent 
pressure in the ID models of stellar atmospheres and stellar evo- 

it appears 



lution, therefore, still remains open. In view of Fig. 15 



unlikely that this can be addressed within the framework of MLT 
when assuming a sharp boundary of the convectively unstable 
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region. In one way or another some "softening" by overshoot- 
ing needs to be introduced. 3D models can give some quantita- 
tive guidance here. The indirect feedback of horizontal inhomo- 
geneities on the vertical structure is also included in this way. 

Obviously, the 3D hydrodynamical stellar atmosphere model 
studied in this work is rather an exploratory one. There is still 
some way to go for improving the 3D models so that they could 
match the realism of the current ID stationary stellar model at- 
mospheres in the treatment of opacities, equation-of-state, ra- 
diative transfer, and scattering. We believe, nevertheless, that the 
advantages of the 3D hydrodynamical models over their classi- 
cal ID counterparts are obvious, which warrants further efforts 
towards their firmer implementation in the field of stellar atmo- 
sphere studies. 
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Appendix A: Computation of the characteristic time 
scales 



The time scales in the red giant explored in this work (see Sect. |3.2.3) were 
calculated Following the prescriptions given in Ludwig et al. ( 2002 2006). 

The radiative relaxation time, £ rac i, can be evaluated using the MLT Formula 



^rad 



PCpAr e 



H)- 



(Al) 



where p is gas density, c p the specific heat at constant pressure, A = q^mlt^p 
the mixing-length, <x SteFan-Boltzmann constant, T gas temperature, f$ = 16 
and fy = 2 dimensionless constants in the MLT Formulation oF Mihalas (see 
Ludwig et al. 1999). The optical thickness oF a convective element, r e , is defined 
as r e =^pA, where x is opacity. 

It can be shown that in the c ase o F non-grey radiative transFer with opacities 
grouped into i opacity bins, Eq. jAlJ can be re- written as 



^rad : 



pc p A 

/30-T3 



/4 + rfJ 



(A.2) 



where the optical thickness oF convective element is now difFerent in each opacity 
bin 



Ti =*ipA 



(A.3) 



and^i is opacity in bin i. Eq. JA2J was used to estimate the radiative time scale 
in the atmosphere oF a red giant model studied here, with the weights w\ For the 
difFerent opacity bins evaluated using 



T db { 
h+ 4df 



(A.4) 



where b[ = BJB and B = ^T A . The mixing-length parameter used in the calcu- 
lations was set to c*mlt=1-8, in accordance with the value derived in Sect. |3.2.5] 
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The advection time scale, £ a dv, was estimated as time interval during which 
the convective element travels one mixing-length 

'adv = -=;, (A.5) 
v z 

where v™ s is temporally and horizontally averaged vertical RMS velocity. 
The adiabatic Brunnt-Vaisala period, fey, was obtained using 



with 



In 



4v= ^"(Vad-V). 



(A.6) 



(A.7) 



where 8 = -(ff )p is the thermal expansion coefficient at constant pressure, g 
gravitational acceleration, V a d and V the adiabatic and actual temperature gradi- 
ents, respectively. 

The Kelvin-Helmholtz time scale, fan, was calculated using 



(A.8) 



where P is gas pressure. 

Finally, the relative importance of radiative and advective heat transport was 
estimated using the Peclet number 



Pe = 



^rad 
^adv 



Appendix B: Changes of the radiative equilibrium 
temperature due to spatial brightness 
fluctuations 



In Sect. |3.2.2] we found a slight (~ 20 K) increase of the mean atmospheric tem- 
perature of the 3D model in the upper atmosphere relative to the radiative equi- 
librium temperature of a plane-parallel ID model atmosphere. 3D and ID model 
share the same effective temperature (and gravity). We want to demonstrate here 
that this temperature increase in the 3D model can be interpreted as a change 
of the radiative equilibrium temperature caused by the spatial fluctuations of the 
brightness in deeper atmospheric layers due to the convection pattern present in 
the 3D model. To this end, we develop a simple model describing the establish- 
ment of the radiative equilibrium temperature when brightness fluctuations are 
present. 

The radiative equilibrium temperature is established by the balance of ab- 
sorbed radiation coming from deeper layers and the radiation emitted accord- 
ing to the local temperature. We consider the radiation field as simply given by 
the Kirchhoff-Planck function B A (T) of radiation temperature T. In the plane- 
parallel case the balance between emission at local temperature T\ and absorp- 
tion of radiation coming from the deeper atmospheric layers with radiation tem- 
perature Tq can be written as 



dAK A (Ti)B A (Ti) = fa dAK A (Ti)B A {W 



(B.l) 



where k a is the opacity and fa the solid angle subtended by the radiating deeper 
layers. We envision this to be the "stellar surface", and To is close to the ef- 
fective temperature. For brevity, we did not explicitely note the pressure depen- 
dence of the opacity k a (P\,T\) in the above equation. All wavelength integrals 
in this section should be taken from zero to infinity, fa is \ in the case of no 
limb-darkening, decreasing to smaller values as fa = |(1 - jd) for a limb- 
darkening being described by a linear limb-darkening law with coefficient a. We 
now introduce brightness fluctuations by considering a two component model 
for the granulation having equal surface area fractions and radiation tempera- 
tures To + ATo and To - ATo. The total emitted flux should correspond to the 
unperturbed situation so that we have to fulfill the normalization condition 



fB JdA^ 



fB \dX- {B A (T + AT ) + B A (T - AT )} 



dAB A (T ), 



(B.2) 



introducing the normalization factor The temperature fluctuations lead to a 
"hardening" of the radiation field, i.e., a shift of the flux towards the blue part 
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Fig. B.l. Solid line: change of the radiative equilibrium temper- 
ature AT i due to horizontal temperature fluctuations according 
Eqs. ( |B.2| ) and ( |B.4| ). Dash- dotted line: dependence as given by 
the linearized solution (IB. 5b. Further details see text. 



of the spec trum. Since JdAB A (T) = ^T A (cr indicates Stefan-Boltzmann's con- 
stant), Eq. (B.2) could be solved exactly for ATo. However, here we are content 
with an approximate - and much simpler - solution for the case ATo ^ To and 
obtain to leading order 



(A.9) fy; 



1 



(B.3) 



We now evaluate the change in the radiative equilibrium temperature AT\ 
relative to the plane-parallel situation due to brightness f luctu ations associated 
with temperature changes of size ATo. We use the ansatz (B.l) to obtain for the 
relation between emission and absorption at a location in the optically thin layers 
of the atmosphere 



dAK A B A (T 1 +AT 1 ) = f B fa 



c1Ak a ^{B a (To + ATo) + B a (T - 



Expanding for AT\ <tc T\ results to leading order in 



AT\ 



fdAK^) 


r 1 d 2 B A 
12 dTl 


T 

To 


I ~ B A (T )\ 


fdA k a (T\) 


dB A 
dT 





ATq)}.(B.4) 



(B.5) 



Equation (B.5) shows that the influence of horizontal inhomogeneities on the ra- 
diative equilibrium temperature is small since AT\ goes as the square of the tem- 
perature fluctuations. Moreover, AT\ is controlled by the complex dependence 
of the opacity k a on wavelength: AT\ vanishes if the opacity is wavelength- 
independent. To ac count for the wavelength dependence of the opacity the in- 
tegrals in Eq. (B.5) have to be evaluated numerically. In do ing s o, w e deci ded 
to go one step back and evaluate the fundamental equations (B.2) and (B.4) be- 
fore expanding for small ATo and AT\ . Moreover, we used the binned opacities 
as applied in the radiative transfer of the 3D and ID model to approximate the 
integrals improving consistency of our approximate analytical model with the 
detailed models. A solution for AT\ for given ATo can be found iteratively. For 
completeness, we also calculated the dependence according Eq. (B.5) , again, 
evaluating t he in tegrals using the binned opacities. 

Figure \BA\ depicts the result for T = 3660 K, T\ = 2900 K, and Pi = 
lOOdyn cm . To corresponds to the effective temperature of the detailed nu- 
merical models, T\ and Pi to the conditions prevailing around optical depth 

logTRoss 3. T\ was adjusted by setting the solid angle fraction to fa = 0.28 

(indicating a rather strong limb-darkening which is actually present). The model 
indeed predicts a heating of the outer atmosphere while a cooling would be also 
possible depending on the wavelength-dependence of the opacity. A white light 
intensity contrast of about 18 % as suggested by Fig. [^corresponds to a tempera- 
ture fluctuation of 4.5 %. At this level of fluctuations our simple analytical model 
predicts an increase of the radiative equilibrium temperature in the optically thin 
layers close to the ~20K found in the 3D model (see Fig.[T3). Hence, we con- 
sider it plausible that the the temperature increase in the optically thin layers in 
the 3D model relative to ID plane-parallel models is a consequence of horizontal 
brightness fluctuations associated with granulation. 
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